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ABSTRACT 

Column density distribution functions of C IV with 12.05 < log(A^) < 14.35 
and Si IV with 11.70 < log(A^) < 13.93 have been obtained using 81 C IV 
absorbers and 35 Si IV absorbers redward of the Lya forest in the lines of sight 
to seven quasars with 2.518 < < 3.78. These distribution functions have 
been directly integrated to yield ion densities at z = 3 to 3.5 of Qq iv = (2-0 ± 
0.5) X 10"^ and Qsi iv = (7.0 ± 2.6) x 10^^ with Hq = 65 km s"^ Mpc"^ and 
go = 0.02 (1 a errors). A larger sample of 11 quasar lines of sight was used to 
measure C II/C IV, Si Ill/Si IV, and N V/C IV ratios, which suggest that C IV 
and Si IV are the dominant ionization stages and that corrections to ficarbon and 
^Silicon are no more than a factor of two. Normalizing the a-process elements 
to silicon and the Fe-copro duct ion elements to carbon gives a density of heavy 
elements in these forest clouds of f^metais = (3.3±0.8) x 10"'' (Hq = 65, go = 0.02). 
The implications for the amount of star formation and for the ionization of the 
IGM prior to z = 3 are discussed. 



Subject headings: early universe — galaxies: formation — intergalactic medium 
— quasars: absorption lines 
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1. Introduction 

In any CDM or CDM + A cosmology, the first generation of stars is expected to form at 
redshifts 10-30 in small (~ 10^ Mq) bound objects and is assumed to imprint the intergalactic 
medium by both enriching and ionizing it. This enrichment may now have been detected in 
Keck spectra of the Lyman forest (Cowie et al. 1995; Tytler et al. 1995; Songaila & Cowie 
1996 [SC96]), and various attempts have been made to use the first crude measurements of the 
forest metallicity (SC96; Haehnelt, Steinmetz, & Ranch 1996; Ranch, Haehnelt & Steinmetz 
1997) to place upper limits on the amount of early star formation in small galaxies (SC96; 
Miralda-Escude & Rees 1997; Gnedin & Ostriker 1997; Haiman & Loeb 1997) and to infer 
the rate of the high redshift type-II supernovae that may be our best hope of detecting 
objects in these early stages of the galaxy formation era (Miralda-Escude & Rees 1997). 

The critical quantity for normalizing the amount of early star formation is the density of 
metals at early epochs. This has generally been obtained (e.g. Miralda-Escude & Rees 1997; 
Haiman & Loeb 1997) by combining the rough estimates of metallicity in forest clouds with 
the total baryon density {fib) inferred from big bang nucleosynthesis (SBBN). Apart from the 
uncertainty in the value of fib from SBBN, this methodology assumes a uniform enrichment of 
metals in the IGM and also entails considerable uncertainty in the determination of individual 
cloud metallicities caused by poorly known ionization corrections and abundance patterns 
in the heavy elements. 

It is possible to avoid many of these problems by directly integrating the observed ion 
column densities in the forest clouds to obtain flion and, within plausible limits, fimetais- 
Carbon and silicon are both well suited to this approach: as discussed below, their dominant 
ionization stages in the IGM are C IV and Si IV, strong doublets that are easy to find and 
measure with 10 m class telescopes (§0); they both have strong accessible lines in other 
ionization stages with which to assess ionization corrections (^); and the ratio of a-process 
silicon to Fe-coproduction carbon is useful for determining the abundance pattern at early 
epochs and to give some idea of the nature of the initial mass function (§§). This Letter 
reports the results of such a direct integration. 

2. Observations 

The data used comprises high resolution observations of eleven quasars with 2.518 < 
Zem < 3.78 made with the HIRES spectrograph on the Keck I telescope, for a variety of 
programs, between 1994 November and 1997 April. All spectra have a resolution R = 36, 000 
and have variable wavelength coverage in the range 3500 A — 7000 A. Full details of the 
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extraction and processing of these spectra are given in a companion paper (Songaila 1997). 
Seven of the spectra have complete wavelength coverage between the quasar's Lya and C IV 
emission lines, whereas the remainder have only partial coverage redward of Ly« emission 
or have remaining inter-order gaps from HIRES's incomplete coverage above 5000 A. 

The sample used to construct the C IV column density distribution consists of all C IV 

doublets detected redward of Lya emission in those spectra with complete coverage. The 
doublets were found by inspection of the spectra, and confirmed by consistency of the col- 
umn density and velocity structure in the two members of the doublet. Systems within 
5000 km s"^ of the quasar's emission redshift were excluded from the sample to avoid in- 
cluding proximate systems which have ion ratios dominated by photoionization from the 
quasars themselves. Also excluded were 10 specifically targeted partial Lyman limit systems 
(PLLS) in 9 quasars that were chosen to be observed in a program to look for suitable sys- 
tems in which to measure the primordial ratio of deuterium to hydrogen. The final sample 
consists of 81 C IV doublets, with 2.02 < z < 3.54 and 12.05 < log(A^c iv) < 14.35, free 
from proximity or observational selection bias. 

The Si IV column density distribution was similarly determined from all Si IV doublets 
detected in the same subsample of quasar spectra as used for the C IV distribution. Only 
systems blueward of the quasars' rest-frame Si IV but redward of the Lyman forest were 
included, again excluding proximate systems and targeted PLLS. This resulted in 35 systems 
in seven lines of sight between z = 2.16 and z = 3.54, with 11.70 < log(A^si iv) ^ 13.93. 

In all cases, ion column densities were determined by fitting up to ten Voigt profiles to 
each redshift system, defined in this context to be all absorption near a given fiducial redshift 
with gaps in velocity space of no more than 50 km s~^. The column density at a given redshift 
is then the total column density of all such components. In all but a few systems, individual 
lines were unsaturated, so fitted column densities are insensitive to 6-value. 

With these samples, the column density distribution function f{N) was determined 
for each ion. f{N) is defined as the number of absorbing systems per unit redshift path 
per unit column density, where at a given redshift z the redshift path, X{z), is defined as 
X{z) = I [(1 + zf - 1] for go = or by X{z) = | [(1 + zf/^ - 1] for qo = 0.5 (Tytler 1982). 
The C IV and Si IV distribution functions are shown in Figure 1, over the redshift ranges, 
2.02 <^<3and3<^< 3.54 for C IV, and 2.16 <z<3and3<z< 3.54 for Si IV. The 
data were calculated with qq — and are plotted with 1 a error bars calculated from the 
Poisson errors based on the number of systems in each bin. 

Because of the somewhat heterogeneous nature of the parent observations, the signal- 
to-noise of the final sample is quite variable, in some cases because the low quasar emission 
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redshift entailed observing below 4000 A where the CCD efficiency dropoff is quite severe, 
in others because of variable coverage at the red end of the spectrum to fill in inter-order 
gaps. Exposure times were also variable for the usual observational reasons. The amount 
of incompleteness at low column density was assessed by recalculating the column density 
distributions using only systems drawn from the highest signal-to-noise spectrum, toward 
Q1422+231, which has an exposure time of 580 minutes and reaches a 1 a limiting column 
density of 1.3 x 10^^ cm~^ for C IV and 5 x 10^° cm~^ for Si IV for 6 = 6 km s~^. Comparing 
this with the distributions obtained from the full sample, it is found that the turnover at low 
column density {Nq iv < 2.8 x 10^^ cm^^) in the full sample is almost entirely a result of 
incompleteness. To take this into account, power laws were fitted only above 5.6 x 10^^ cm~^ 
for C IV and 10^^ cm~^ for Si IV , with best-fit indices of —1.5 for both low and high redshift 
C IV, and —1.8 for low redshift Si IV and —2.0 for high redshift Si IV. These values are quite 
similar to those measured in H I in the forest at these redshifts (Petitjean et al. 1993; Hu et 
al. 1995; Kim et al. 1997) and in higher column density C IV samples (Petitjean & Bergeron 
1994). 

Qc IV and Qsi iv were calculated from the column density distributions of Figure 1 
according to 

O _ Hp X] ^ion /-.N 
''ion A~\7 ^ion {^J 

Pcrit AA 

where pcrit = 1-89 x lO^^^/i^ g cm"^ is the cosmological closure density, mion is the ion's 
mass, and Hq = 100 h km s^^ Mpc^^ (e.g., Lanzetta et al. 1991). Values of fic iv and 
^Si IV were calulated for z > 3 in individual lines of sight to the five quasars in the sample 
with > 3 and complete spectral coverage between the quasar's Lya and C IV emission 
lines. The results are tabulated in Table |l| for Hq = 65 km s~^ Mpc^^ and qo = 0.02 
(f2ion scales as [1 + 2qozY^'^) and give some idea of the uncertainty in calculating ilion- 
Formal mean values, weighted by AX, are flc iv = (1.2 ± 0.3) x 10~^ [1 + 2qozY^'^ 
and ilgi IV = (4.3 it 1.6) x 10~^ h"^ [1 + 2go^]^^^ , where the errors are 1 a. The mean 
redshift of the sample is (2;) = 3.18. A similar procedure applied to the z < 3 systems gives 
ric IV = 9.3 X 10-^ h-^ [1 + 2qozY/^ and figi iv = 5.0 x 10"^ h'^ [1 + 2go^]^/^ at {z) = 2.57. 

Because of the rarity of high column density systems {Nc iv ^ 3 x 10^^ cm^^) consid- 
erably longer path lengths are required to determine the number density of such systems, 
and the present data measure f2ion only for absorption with column density less than this 
value. (This corresponds roughly to Nn i < 10^'' cm~^, from SC96.) The metal densities of 
the stronger systems are addressed in a companion paper (Songaila 1997). The contribution 
of systems to Qq iv and iv converges at the low column density end, and systems weaker 
than those observed will not contribute significantly to the density unless there is a very 
rapid upturn below the observed range. 
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3. Ionization Balance 

A less restricted sample was used to assess the ionization balance correction to be applied 
to the C IV and Si IV sample. This was drawn from the full set of eleven quasar spectra, 
excluding proximate systems. In cases where there was no apparent absorption at the C IV 
redshift in the other ionic species, an upper limit to the absorption in other ions was found 
by formally fitting the C IV component model to the local continuum. This approach results 
in conservative upper limits, especially to lower ionization species, since C IV absorption is 
likely to be more widespread in velocity space than lower ionization absorption. 

As is shown in Figure 2, C II is a trace ion relative to C IV with nearly all systems 
having C II ^ 0.1 C IV. An additional direct search in the spectra for C II systems with no 
strong C IV yielded no additional systems. For photoionization models, this in turn yields 
a high photoionization parameter (logF ^ —2) which, for a wide range of photoionizing 
spectra, implies that Si III ~ Si IV and C III ~ C IV (e.g., Steidel 1990; Giroux & Shull 
1997). For Si III this can be directly verified from observations (middle panel of Figure 2) 
which show that, even with Lja forest contamination, the upper limits on the Si III column 
density are generally compatible with Si IV measurements. Ionization by a starburst galaxy 
spectrum with no high energy photons could result in a much higher ratio of C III to C IV 
while maintaining the low C Il/C IV and Si Ill/Si IV ratios, but it would also produce Si IV 
much in excess of C IV, which is not observed: at 2; > 3 Si IV/C IV has an average value of 
about 0.2 (SC96; Songaila 1997). 

Finally, neither N V nor O VI is strong in the C IV-selected systems. The right panel of 
Figure 2 shows that N V < 0.05 C IV in the small number of systems measured. The high 
observed ratio of Si IV to C IV combined with the low observed values of C Il/C IV and 

N V/C IV suggests that forest clouds at 2; > 3 are ionized by a broken power low spectrum 
with relatively few high energy photons and that it is unlikely that there is much material 
above the C IV and Si IV levels (SC96; Giroux & Shull 1997). The best guess, therefore, 
would give l^carbon 2 IV a^d f^siiicon ^ 2 ly in thcsc clouds. 

4. Discussion 

Converting ^carbon and ^silicon to the fimetais of all metals requires an assumption of 
an abundance pattern in the z > 3 forest. As in the old metal-poor halo stars, sili- 
con is overabundant with respect to carbon relative to solar. However, assuming that 
the carbon and silicon abundances trace the universal abundances of iron coproduction 
and a- process elements, respectively (e.g., Timmes, Lauroesch, & Truran 1995), then for 
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the Fe-coproduction elements (C, N, Fe), Q — 1.75 x Qc, and for the a-process elements 
(O, Ne, Si, Mg, S), = 18.5 x f^si, and setting Qc = iv and Qsi = 2^31 iv gives 
^metals = (2.1 ± 0.6) X 10~'^ [1 + 2go-2]^^^ (relative abundances from Anders & Grevesse 
1989). Assuming a value of from SBBN of 0.005 < n^h^ < 0.016 (Songaila, Wampler, 
& Cowie 1997) gives l^mctais/^^b = (1-3 ± 0.4 ^ 4.2 ± 1.2) x 10"^ [1 + 2qozY/^. With 
Hq = 65 km s~^ Mpc^^, Qq = 0.02, and a solar metallicity of 0.019, this implies a minimum 
universal metallicity relative to solar, at 2; ~ 3, in the range [—3.3] to [—2.8]. (The value of 
D/H ^ 2.3 X 10-^ ± 3 X 10^^ (statistical) ± 3 x 10"^ (systematic) given in Tytler, Fan & 
Buries (1996) has now been revised upwards to 3.3 x 10~^ (Tytler 1997) which would imply 
Jlfe h'^ — 0.016 and a metallicity of [—3.3], while the value of 2 x 10~^ obtained by Webb et al. 
(1997) would give a metallicity of [—2.9] [Hq = 65, go = 0.02].) This is a minimum range of 
metallicity since there may be additional metals in higher column density clouds or galaxies 
as well as in ionization states that were not sampled by these observations. Furthermore, 
adoption of go = 0.5 would roughly double the metallicity. However, much of Qi, inferred 
from SBBN is believed to reside in the forest clouds at this time (Ranch et al. 1997; Kim 
et al. 1997), suggesting that this calculated metal density should be a good estimate of the 
total metals. Interestingly, since a metalhcity of 10~^ gives rise to one ionizing photon per 
baryon (e.g., Miralda-Escude & Rees 1997), the metallicity measured here is just sufficient 
to preionize the IGM. However, this metallicity range is substantially lower than the value 
of 2.4 X 10^"^ assumed by Miralda-Escude & Rees (1997) and Haiman & Loeb (1997) and 
implies that their predicted rate of very high-;^ type-II supernovae could be lowered by as 
much as an order of magnitude, to a value of around 1 SN per 10 arcmin^ per year. 

I am grateful to the many people at the Keck telescopes who made these observations 
possible, and to Len Cowie and Esther Hu for obtaining some of the observations on which 
this work is based. The research was supported by the National Science Foundation under 
grant AST 96-17216. 



- 7- 



REFERENCES 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Cowie, L. L., Songaila, A., Kim, T.-S., & Hu, E. M. 1995, AJ, 109, 1522 
Giroux, M. L., & Shull, J. M. 1997, AJ, 113, 1505 

Gnedin, N. Y., & Ostriker, J. P. 1997, ApJ, 486, in press [preprint |astro-pii/9612127[ 

Haiman, Z., & Loeb, A. 1997, in Proceedings of the 18th Texas Symposium on Relativistic 
Astrophysics, ed. A. Ohnto, J. Frieman, & D. Schramm (Singapore: World Scientific), 
in press [preprint |astro-ph /9 70 1 239[| 

Haehnelt, M. G., Steinmetz, M., & Ranch, M. 1996, ApJ, 465, L95 

Hu, E. M., Kim, T.-S., Cowie, L. L., Songaila, A., & Ranch, M. 1995, AJ, 110, 1526 

Kim, T.-S., Hu, E. M., Cowie, L. L. & Songaila, A. 1997, AJ, 114, 1 

Lanzetta, K. M., McMahon, R. G., Wolfe, A. M., Turnshek, D. A., Hazard, C, & Lu, L. 
1991, ApJS, 77, 1 

Miralda-Escude, J., & Rees, M. J. 1997, ApJ, 478, L57 

Petitjean, P., Webb, J., K., Ranch, M., Carswell, R. F., & Lanzetta, K. 1993, MNRAS, 262, 
499 

Petitjean, P., & Bergeron, J. 1994, A&A, 283, 759 

Ranch, M., Miralda-Escude, J., Sargent, W. L. W., Barlow, T. A., Weinberg, D. H., Hern- 
quist, L., Katz, N., Gen, R., & Ostriker, J. P. 1997, ApJ, in press [preprint |astro- 
ph/96l^^ 

Ranch, M., Haehnelt, M. G. & Steinmetz, M. 1997, ApJ, 481, 601 

Songaila, A. 1997, in preparation 

Songaila, A., & Gowie, L. L. 1996, AJ, 112, 335 [SG96] 

Songaila, A., Wampler, E. J., & Gowie, L. L. 1997, Nature, 385, 137 

Steidel, G. G. 1990, ApJS, 74, 37 

Timmes, F. X., Lauroesch, J. T., & Truran, J. W. 1995, ApJ, 451, 468 
Tytler, D. 1982, Nature, 298, 427 

Tytler, D., Fan, X.-M., Buries, S., Gottrell, L., Davis, G., Kirkman, D., & Zuo, L. 1995, in 
QSO Absorption Lines, ed. G. Meylan (Heidelberg: Springer), 289 

Tytler, D., Fan, X.-M. & Buries, S. 1996, Nature, 381, 207 



-8- 



Tytler, D. 1997, talk presented at the 13th lAP Astrophysics Colloquium, Structure and 
Evolution of the Intergalactic Medium from QSO Absorption Line Systems, Paris, 
1-5 July 1997 

Webb, J. K., Carswell, R. F., Lanzetta, K. M., Ferlet, R., Lemoine, M., Vidal-Madjar, A. & 
Bowen, D. V. 1997, Nature, 388, 250 



This preprint was prepared with the AAS IM^jX macros v4.0. 



-9- 



Tablc 1. z > 3 Hioii by Quasar 



Quasar 




AX 


EN{CW)/AX 


EA^(SiIV)/AX 


^CIV 


^SilV 








(xlO^^ cm-2) 


(xlO^^ cni-2) 


(xlO-8) 


(xlO-9) 


0014+813 


3.384 


1.29 


0.62 


0.42 


1.2 


1.8 


0302-003 


3.286 


0.88 


0.51 


0.17 


0.9 


0.74 


0956+122 


3.301 


0.94 


2.4 


3.7 


4.5 


16.0 


1159+123 


3.493 


1.76 


0.62 


0.63 


1.2 


2.7 


1422+231 


3.620 


2.32 


1.3 


2.7 


2.4 


12.0 


Mean 






(1.1 + 0.3) 


(1.6 + 0.6) 


(2.0 + 0.5) 


(7.0 + 2.6) 



Note. — For all five quasars, Si IV and C IV lie above the Lya forest at z > 3 so AX 
is the same for C IV and Si IV. Q is computed for Hq = 65 km s~^ Mpc~^ and Qq — 0.02. 
Mean values are weighted by AX and quoted errors are 1 a. 
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Log N(C IV) Log N(C IV) 

Fig. 1(a). — C IV column density distribution functions for 2; < 3 and 2; > 3 for a total of 
81 absorption line systems in seven lines of sight toward distant quasars. The bin size is 
]\[ cm^^. Error bars are ±1 a, based on the number of systems in each bin. The solid 
lines show the best fit power law with index a = —1.5, fitted to data above 6 x 10^^ cm~^ 
where the data is substantially complete. 
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Log N(Si IV) Log N(Si IV) 

Fig. 1(b). — As in (a) for Si IV, with a total of 35 systems and a bin size of 10'''^ cm~^ 
for z < 3 and 10^''^ N cm^^ for z > 3. The soUd hnes are power laws of index a — —1.8 for 
z < 3 and o; = — 2 for 2; > 3, fitted to data between 10^^ cm~^ and 10^^ cm~^. 
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Fig. 2. — Ion ratios for absorption line systems toward the full quasar sample. The C II and 
N V samples were C IV-selected (A^c iv > 1-5 x 10^^ cm'^ for C II; A^c iv > 4 x 10^^ ^.^-2 
for N V) and the Si III sample was Si IV-selected (A^si iv > 2 x 10^^ cm-^). For C II and 
N V, the sample was restricted to lines lying longward of the Lya forest. This is not possible 
for Si III A1206 A, where all measured values were treated as upper limits because of forest 
contamination. 



